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The numerical simulations of an ultrashort pulse propagation in a one-dimensional nonlinear 
photonic crystal are carried out. It is known that the relaxation of cubic nonlinearity is the reason 
for the effect of pulse self-trapping in such multilayer system. In this paper we study further 
implications of this effect. It is shown that the trapped light absorbs additional low-intensity pulses 
which cannot be self-trapped per sc. On the other hand, such low-intensity pulses are subject of 
the so-called induced trapping when light gets trapped due to a collision of two such pulses. We 
consider the conditions for this effect in both cases of co- and counter-propagating pulses. 
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I. INTRODUCTION 

It is known that nonlinear response of optical media 
is not generally instantaneous and is described by a cer- 
tain settling time. This relaxation of nonlinearity can 
be neglected if the characteristic time of electromagnetic 
field (such as a pulse duration) is much greater than the 
relaxation time. However, in the modern era of ultra- 
short pulses, there is growing number of situations when 
this neglect cannot be justified. The standard means 
to calculate the relaxation of non-resonant cubic (Kerr) 
nonlinearity is the so-called Debye model [l|, 



, d,Sn 

tnl—7- +071 = n2l, 

at 



(1) 



where Sn is the nonlinear part of refractive index, n2 the 
Kerr nonlinear coefficient, and tni the relaxation time. 
The latter depends on the specific mechanism of nonlin- 
earity. In this paper it is assumed to be of the order of 
several femtoseconds. 

The influence of non-instantaneousness of nonlinearity 
on optical response was studied in a number of works dur- 
ing last several decades. Not aiming to name all of them, 
we can call researches on laser beam self- focusing [I-HIj 
filament formation [F] , parametric amplification [bf , pulse 
compression. [7] , modulation instability .Sr.llJ , pulse train 
generation [I^], soliton-array generation [IJl, instab ility 
of speckle patterns fl^ , so litary pulse dynamics [H, [l6| . 
optical switching llSl j. etc. In this paper we deal 
with ultrashort (femtosecond) pulse propagation in a one- 
dimensional nonlinear photonic crystal which is a set of 
periodically arranged dielectric layers. Such multilayer 
structure can be symbolically designated as (AB)^ {N 
is the number of pairs of layers of A and B type). One 
of the first studies of role of nonlinearity relaxation in 
such a system is the paper by Vlasov and Smirnov fl^ 



where the problem of pulse com pres sion was under in- 
vestigation. In recent works |20l - [23| the effect of pulse 
self-trapping in the photonic crystal due to relaxing non- 
linearity was discovered and discussed. The present pa- 
per is a continuation of those researches, so that it is 
worth recalling briefly the main results obtained there. 

It was shown that the pulse of high enough inten- 
sity can be trapped inside the one-dimensional photonic 
crystal due to formation of dynamical nonlinear " cavity" , 
or " trap" . This trap appears only if both linear refrac- 
tive index modulation and Kerr nonlinearity relaxation 
present. The range of pulse durations and relaxation 
times for such self-trapping effect to be observed in the 
multilayer system of several hundreds layers was stud- 
ied as well: i„; varies from a fraction of a femtosecond 
to more than 100 fs, and tp from about 10 fs to sev- 
eral hundreds fs. This corresponds to the fast electronic 
mechanism of Kerr nonlinearity. Though we do not mean 
any specific materials, it is believed that such nonlinear 
structures can be composed, for example, from doped 
glasses with rapidly relaxing nonlinearities. This leads 
to the requirement of comparatively high intensities of 
the pulses because of the well-known from experiments 
approximate proportionality 712 ^ tni [H. Nevertheless, 
the necessary peak intensities of the order of a hundred 
of GW/cm^ seems to be not excessively high for fem- 
tosecond pulses from the viewpoint of optical damage of 
the materials. Polymeric materials are worth to be men- 
tioned due to both high nonlinearity coefhcient and fast 
relaxation (23 | that make them perspective for applica- 
tions in the discussed situations. In [2^ we studied in 
detail the conditions for self-trapping in different con- 
figurations of the structure (linear and nonlinear layers, 
focusing and defocusing nonlinearities, etc.) taking into 
account the correlation between the nonlinearity coef- 
ficient and the relaxation time mentioned above. An- 
other problem studied is the spectral transformations of 
light pulses interacting with nonlinear photonic crystal 
in the regime of self-trapping [2l[. In particular, un- 



' |dvnovitsky@tut.by"] 



der properly chosen conditions, it is possible to gener- 
ate quasi-monochromatic radiation or quasi-continuum 
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covering the whole band gap. Finahy, the possibihty of 
asymmetric light transmission due to pulse self-trapping 
was analyzed recently [23l |. 

In this paper we consider the situation of not a sin- 
gle pulse, but of many pulses inside the photonic crystal 
with relaxing nonlinearity. We are especially interested 
in investigation of interaction of probe pulses with the 
previously trapped radiation. This problem is studied in 
im Another question is connected with possibility of light 
trapping due to a collision of two low-intensity pulses 
which are not trapped when propagate separately. This 
situation, which we call the induced trapping, is consid- 
ered in mil 

II. TRAPPING OF PROBE PULSES 

At first, let us state the main equations used in this 
paper. Propagation of an ultrashort pulse in a one- 
dimensional nonlinear photonic crystal is described by 
the Maxwell wave equation 

where E is the electric field strength, n is the refractive 
index which depends on light intensity / = \E\'^ as 

n{z,t) — nQ{z) + Sn{I ,t). (3) 

Here no{z) is the linear part of the refractive index vary- 
ing periodically along the z axis, and the nonlinear ad- 
dition 5n{I,t) behaves according to Eq. ([1]). Further 
we consider femtosecond light pulses of Gaussian shape 
with the amplitude A = Ajnexp{—t'^/2tp), where tp is 
the pulse duration. To analyze the interaction of such 
a pulse with a nonlinear photonic crystal, i.e. to solve 
self-consistently the system (HJ-dS]), we use the finite- 
difference time-domain method of numerical simulations 
which was described in detail in Ref. ^2^. The stability 
of the algorithm used is governed by the well-known von 
Neumann condition At/Az < 1/v ~ n/c, where At and 
Az are the time and space steps, respectively. If the ratio 
of the steps is small enough (taking into account possi- 
ble nonlinear change of n), one provides the necessary 
stability of calculation. The level of discretization (size 
of the steps) is chosen to be optimal in respect to both 
accuracy and calculation time and allows to obtain the 
reliable general dependencies discussed further. 

The parameters of the photonic crystal [the structure 
of (AB)^ type] used in our calculations are as follows: 
the linear parts of refractive indices of the layers A and B 
are Ua ~ 2 and ni, — 1.5, respectively; their thicknesses 
a — 0.4 and b — 0.24 fim; the number of layers N — 
200. The pulse duration is tp — 30 fs, and the central 
wavelength of the initial pulse spectrum is Ac — 1.064 
fim if not stated otherwise. The nonlinear coefficient of 
the material is defined through the nonlinear term of the 
refractive index, so that 712/0 = 0.005; this means that 
the pulse amplitude is normalized by the value ^0 = 



^/Tq. The relaxation time of the nonlinearity of both 
layers is tni — 10 fs. We adopt these parameters in this 
paper though the similar effects can be obtained even 
in the half-linear structure [23|. It is also important to 
note that the parameters (in particular, the wavelength) 
satisfy the requirements on the sign of group velocity 
dispersion studied in Refs. [21I. [23|. 

Figure [ij a) shows the dependence of the output energy 
leaving the photonic crystal during a certain time inter- 
val after pulse incidence on the pulse amplitude. The 
output energy is calculated by integration of intensity of 
light leaving the photonic crystal over time. The relative 
output energy normalized by the input one is 

^ ^ g-* ^ Ji°-\t)dt 

Qinc J /mc(i)^^' > 

where Q™"^ and Q""* are the absolute values of energy of 
the incident and output light, respectively; /""^ and /°"* 
are the corresponding intensity profiles. The value cal- 
culated by Eq. (|4]) at the input edge (incidence plane) of 
the structure gives the normalized energy of the reflected 
light, while the energy calculated at the exit edge corre- 
sponds to the transmitted light. The sum of these two 
values is the total output energy. The relative output en- 
ergy in Fig. [Ija) is integrated over the time 300ip which 
is approximately ten times larger than the pulse trans- 
mittance time in the linear regime (about 30ip). The dip 
seen in the curve for the total output energy is the feature 
of self-trapping of the pulse. This means that the pulses 
with the amplitudes in the certain range are trapped in- 
side the photonic crystal. Outside this range, the low- 
intensity pulses transmit through the system, while the 
high-intensity ones are mostly reflected. 

The result of self-trapping is the formation of localized 
light intensity distribution which is the indication of trap 
creation. This trap stores most of the pulse energy for 
the times of the order of several thousands tp [20]. It 
is worth noting that self-trapping does not mean light 
absorption, because the modulated refractive indices re- 
main real. Let us study what happens when the second 
(probe) pulse interacts with this excited state of the pho- 
tonic crystal containing trapped light. In this case the 
intensity of light in Eq. [1] is governed by the sum of the 
field present at a certain space point and time instant. 
The second pulse starts at the instant lOOip after the first 
(trapped) one. The results of calculations of the output 
energy (integrated over the time 300ip) as a function of 
amplitude of the probe pulse is presented in Fig. [Ub). 
It is seen that the trap formed due to the pulse with 
the amplitude Am = iA^ blocks propagation of probe 
pulse with low enough intensity. This means that the 
energy of the probe pulse gets stored inside the photonic 
crystal, so that the trap collects more and more light. 
Our estimate shows that the probe pulse with Am = Aq 
loses about 60% of its energy due to interaction with the 
trap. The region of optimal self-trapping of single pulses 
(3^0 < Am < 4:Ao) is naturally the range where the most 
part of both pulses is trapped. Finally, the high-intensity 
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FIG. 1. (Color online) Dependence of the output light energy 
(normalized to the input energy) (a) on the peak amplitude 
of the incident pulse, (b) on the peak amplitude of the probe 
pulse interacting with the photonic crystal in the excited state 
(after trapping of the initial pulse with the amplitude Am — 
3Ao). The energy is integrated over the time SOOtp. 



probe pulse (with Am > 5Ao) cannot be trapped. This, 
however, is connected not with the breakdown of the trap 
but with reflection of the high-intensity pulses seen in 
Fig. ma). 

We prove these conclusions in Fig. [2] where the distri- 
butions of light intensity inside the photonic crystal at 
different time instants are shown. It is seen that, by the 
time t — lOOtp (launch time of the probe pulse), the trap 
(bell-shaped stable light distribution) formed by the first 
pulse with the amplitude Am = 3Aq exists inside the mul- 
tilayer system (the total length of the system is 128 /im). 
If the probe pulse has low peak amplitude [Am = Aq, 
see Fig. Ela)], it is just absorbed by the trap, the peak 
intensity of the distribution increasing from about 12Iq 
to 14/o. This distribution is slowly spreading with time. 
The situation is totally different for the probe pulse with 
the amplitude TAq [Fig. El^b)]. It is seen that this high- 
intensity pulse is stopped near the very entrance of the 



FIG. 2. Distribution of light intensity inside the photonic 
crystal at different time points. The probe pulse amplitude 
is (a) Am. ~ Aq, (b) Am ~ 7Ao. The peak amplitude of the 
first pulse is Am = 3^o- 



photonic crystal and then effectively reflected. However, 
some part of its light penetrates in the vicinity of the 
trap and perturbs it. As a result of this perturbation, 
the peak intensity of the distribution decreases (about 
8/o at 300ip), the trap spreads and shifts towards the 
entrance of the structure. At time t = lOOOtp the distri- 
bution is strongly widen and has maximum of the order 
of 0.5/o. In other words, the trap looses the light energy 
faster than in the case of low-intensity probe pulse and 
is less stable. 

Nevertheless, we can state that, generally, the trap per- 
sists and cannot be overcome by a single probe pulse. To 
study stability of the trap further, we launch more probe 
pulses inside the photonic crystal. The first probe pulse 
starts at t = 50ip after the initial one (which forms the 
trap); the interval between every subsequent probe pulses 
is lOtp. The integral output energy as a function of the 
number of probe pulses is plotted in Fig. Elja). It is 
seen that, at first, the part of energy leaving the mul- 
tilayer structure is growing quite rapidly. But after the 
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FIG. 3. (Color online) (a) Dependence of the output light 
energy (normalized to the input energy) on the number of 
probe pulses. The amplitude of the probe pulses and of the 
initial one is Am = 3Ao; the energy is integrated over the 
time SOOtp. (b) Transmitted intensity profile and (c) intensity 
distribution inside the photonic crystal (at time 300tp) for the 
case of 11 probe pulses. 
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FIG. 4. Distributions of light intensity inside the photonic 
crystal for different number of probe pulses. The peak am- 
plitude of the initial and probe pulses is Am ~ SAq. The 
number of probe pulses and the time points are indicated on 
the panels. 



shift of the trap towards the input of the system. Obvi- 
ously, this process continues as we launch more and more 
pulses in the structure. Such situation is shown in Fig. 
m We also see substantial spreading of the distribution 
for larger number of input pulses. Finally, for 100 probe 
pulses (lower three panels in Fig. the bell-shaped dis- 
tribution tends to collapse sooner than for lesser number: 
the trap is absent at t = 15QQtp which is about 500ip after 
the last pulse enters the system, while for 25 pulses the 
trap is still stable after 700tp (see left upper panel). Nev- 
ertheless, even at i = 2000tp approximately 40% of the 
energy of 100 pulses remains inside the photonic crystal 
due to multiple reflections on the numerous boundaries. 



III. INDUCED TRAPPING 



sixth probe pulse the trapped energy (in relative units) 
stays approximately the same. This means that the ab- 
solute value of light energy trapped inside the system is 
growing. The evidence for this conclusion is presented 
in Fig. [3jc) where the intensity distribution is shown for 
the case of 11 probe pulses. The peak value of intensity 
is about 50/o (at SOOtp) which is much greater than 14/o 
for the single probe pulse [see Fig. [DJa)]. The profile of 
the transmitted radiation shown in Fig. ^h) allows to 
identify the intensity peaks corresponding to particular 
probe pulses. Notice that the first strongly pronounced 
peak is connected with the fourth probe pulse, i.e. the 
first, second, and third probe pulses are almost entirely 
absorbed. Every following pulse appears more or less 
sharply in the output of the photonic crystal. 

Returning to Fig. EJc), one can note not only the in- 
crease in peak intensity of light in the trap but also the 



In Fig. [IJa) we see that the transition between the 
regimes of propagation and trapping is an abrupt one. 
Therefore the amplitude of this transition which is ap- 
proximately equal to 2.2Ao can be called the critical am- 
plitude. The pulses with subcritical intensity freely prop- 
agate through the photonic crystal (the most part of light 
is simply reflected and transmitted), while the supercrit- 
ical pulses are self-trapped. In this section we study the 
possibility for light trapping as a result of collision of 
two pulses with subcritical intensities. We refer to such 
a situation as the induced trapping. 

The first possible scheme is the collision of two co- 
propagating pulses. This scheme implies that the second 
pulse must move faster than the first one and overtake it 
at a certain point of space. As a result, the light intensity 
is to increase in the vicinity of this point and formation 
of the trap is expected due to summation of the fields 
of two subcritical pulses. To demonstrate possibility of 
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FIG. 5. Spectral curves for reflectivity (upper panel) and 
group velocity (lower panel) of the photonic crystal under 
consideration. The parameters of the structure are given in 
the text. 



such course of events, we should select the proper values 
of the parameters of the pulses, so that they could meet 
inside the photonic crystal. It is known that the group 
velocity Vgr of the pulse in the periodic structure depends 
on the wavelength: it is large in the transmission band 
and decreases towards the band gap where it vanishes. 
The reflectivity spectrum and group velocity behavior 
near the band edge of the system considered in this paper 
are shown in Fig. [SJ It is seen that the group velocity at 
the central wavelength Ac = 1.064 fim is about half of the 
light speed {vgr = 0.5c). Therefore, we keep the central 
wavelength of the second ("fast") pulse unchanged (Ac2 = 
1.064 fim) and vary the central wavelength Ad of the first 
("slow") pulse. 

Note that in our consideration, we neglect the disper- 
sion of the materials of the layers, i.e. the change of the 
refractive index with frequency. This is justified by the 
very slow rate of this change for many optical materi- 
als (such as glasses) in the s pec tral range of our interest, 
namely near infrared region |25l . [26j . In addition, we take 
only narrow section of this range from about 1.05 to 1.12 
/zm. Thus, the only source of dispersion in the further 
study is the structural dispersion caused by the order of 
the photonic crystal layers. 

The results of calculations of the output energy as a 
function of Ad are presented in Fig. M^a). The time 
interval between the peaks of co-propagating pulses is 
At = 5tp] their amplitudes are equal and sub-critical 
{Ami = = 2^o)- As it was expected, the effective 
trapping occurs for Ad > Ac2 (i-c. when Vgri < Vgr2), so 
that it can be unambiguously interpreted as a result of 
pulse collision. The details of pulse interaction are clar- 
ified in Figs, [ll^b) andlSJc) where the resulting intensity 
profiles and distributions (at 300tp) are plotted. It is seen 
that, when the first pulse moves faster (Ad = 1.05 /im) 
than the second one, there are two sharp peaks at the out- 



put. For the pulses with equal wavelengths (Ad = 1.064 
/im), these peaks are much less pronounced and about a 
half of the energy of the pulses stays inside the system. 
There is the bell-shaped light distribution near the very 
exit of the structure [see Fig. [lljc), upper right panel] 
which indicates the trap formed in this case. This trap 
is perhaps due to multiple reflections from the interfaces 
between layers as well as the reflection from the output 
boundary of the photonic crystal. This is corroborated 
by the position of the trap (near the very output) and 
the structure of the transmitted radiation (two peaks of 
lowered intensity, in contrast to the single peak seen in 
the case of optimal induced trapping at Ad = 1.075 /xm). 
It is worth noting that the interaction between the pulses 
can be observed when they move close enough one to an- 
other, i.e. the interval between pulses is not very large. 

At Ad = 1.075 /im we have nearly optimal induced 
trapping when there is the single low-intensity pulse at 
the exit. This transmitted pulse corresponds to the first 
incident one, while the second one is almost entirely ab- 
sorbed. The trap forms close to the exit and has com- 
paratively high peak intensity (about 9/o). As we move 
further towards the band gap shown in Fig. [Hfa) by 
the dotted vertical line, the first pulse gets more and 
more reflected, so that the transmittance starts to in- 
crease as well. However, even near the very band edge 
(for Ad = 1.11 /i.m), some part of the light energy (about 
40%, to be exact) is trapped. The corresponding distri- 
bution shows that the trap forms near the entrance of 
the structure. This can be interpreted as a result of the 
collision of the very slowly moving first pulse (or the rest 
part of it which was not reflected) and the second one. 
The resulting distribution has large width and low peak 
intensity (about 2/o). For larger Ad, the output energy 
rapidly increases to the unity as we are now inside the 
forbidden gap. 

Thus, the effect of induced trapping can be observed 
in the rather wide range of wavelengths of the pulse near 
the band edge of the photonic crystal. This range is 
limited, at a given interval between the pulses, by the 
length of the structure and, on the other hand, by the 
reflectance near the band edge. The length of the system 
is important at small differences between speeds of the 
pulses when one of them needs much time to overtake the 
other; at best, they collide near the exit of the system. 
The importance of the band edge appears in the opposite 
case of large velocity difference when, though the pulses 
collide near the entrance of the system, the most of the 
"slow" pulse is reflected. 

In Fig. [7]we study dependence of the induced trapping 
of two identical co-propagating pulses (with amplitudes 
2Aq) on interval At between them. The central wave- 
length of the first pulse (Ad — 1.075 /im) is nearly op- 
timal for trapping. A< = means that, in fact, we have 
a single pulse with the amplitude 4Ao, i.e. the super- 
critical one. The trapping with approximately the same 
efficiency is observed for At < btp but then the output 
energy starts to increase in a stepwise manner. This is 
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FIG. 6. (Color online) (a) Dependence of the output light energy (normalized to the input energy) on the wavelength of the 
first pulse Aci. The amplitude of both pulses is Ami = Am2 ~ 2Ao; the wavelength Ac2 = 1.064 /im; the time-interval between 
the pulses At = 5tp; the energy is integrated over the time SOOtp. (b) Profiles of transmitted and reflected light for different 
Aci. (c) Corresponding distributions of light intensity inside the photonic crystal sd, t — SOOtp. 



accompanied by the rise of transmitted light as the trap 
forms nearer to the exit of the system (see the panels cor- 
responding to At = 7tp and 9tp). Finally, at At = lltp, 
the collision do not happen: the interval is too large for 
the pulses to have time to collide. This is the same limit 
due to the finiteness of the length of the photonic crystal. 

The next question is the dependence of the induced 
trapping on the amplitudes of the pulses at the fixed 
values of other parameters (we adopt At = 5tp and 
Aci — 1.075 /im). Without loss in generality, we can vary 
only the amplitude of the second ("fast") pulse Am2 leav- 
ing the amplitude of the first one the same as previously 
{Ami — 2y4.o). The results of calculations at different 
Am2 are shown in Fig. [5] It is seen in the panel [Sl^a) 
that the output energy starts to decrease at Am2 > Aq 
and reaches the minimum at approximately 2Aq. At 
Am2 > 2.5^0 there are the local peaks in the curves 
of total and transmitted energy. These peaks can be in- 
terpreted as a result of effective trapping of the second. 



high-intensity pulse preventing the collision. This con- 
clusion is proved in Fig. [5fc): the transmitted pulse pro- 
files have approximately the same peak intensity in the 
cases Ajn2 = (the single first pulse) and Am2 = SAq 
(the trapping of the second pulse), while in the inter- 
mediate variant {Am2 = 2Ao) the transmitted pulse is 
strongly suppressed due to the induced trapping. The in- 
tensity distributions in Fig. [H^d) demonstrate the differ- 
ence between the trap formation due to the collision (at 
A,n2 = 2^0 ) and single pulse absorption (at Am2 — SAq). 
Note that, in the last case, the trap forms at larger dis- 
tance than in the case of the single SAq pulse [see Fig. 
[2] at t = 75tp] . This means that the dynamics of the 
second pulse in the two-pulse scheme still depend on the 
first one. 

Figure EJa) does not allow to demarcate unambigu- 
ously the regions of induced trapping and single pulse 
trapping. To carry out this demarcation, we calculate 
the excess in output energy, i.e. the difference between 
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FIG. 7. (Color online) (a) Dependence of the output light energy (normalized to the input energy) on the time-interval between 
the pulses At. The amplitude of both pulses is A,ni = Am2 = 2Ao; the wavelengths Aci = 1.075 ^m, Ac2 ~ 1.064 ^m; the 
energy is integrated over the time SOOtp. (b) Profiles of transmitted and incident light for difi'erent At. (c) Corresponding 
distributions of light intensity inside the photonic crystal at t = SOOtp. 



output energies in the cases of two-pulse and single-pulse 
propagation. This value can be written as 

12 

where the relative output energies W12, Wi, and W2 are 
defined by Eq. (jU in the instances of both pulses, single 
first and single second pulse, respectively. Q™'^, 
and Q™*^ the corresponding absolute energies of in- 
cident light (at that, Q™'^ = Q™= + Qf"). The value 
ISW shows how much extra energy leaves the system 
when both pulses are launched in comparison with the 
single pulse cases; if it is negative, AVK < 0, then one 
can say that the additional energy is trapped inside the 
multilayer structure due to the interaction between the 
pulses. The excess value AMK extracted from the data in 
Fig. El^a) is shown in the panel (b). This figure is the evi- 
dence that, in the region Aq < Am2 < 2.7Ao, the effective 
trapping occurs due to the collision of the pulses (up to 
80% of the total energy can be additionally trapped) . At 



higher values of the amplitude, Am2 > 2.7Ao, the excess 
value is small and approximately constant, AW > —0.1, 
which means that only a minor part of light is trapped 
due to the interpulse interaction. This proves our previ- 
ous conclusion about practically independent trapping of 
the second pulse in the high-intensity regime. 

Finally, we should consider the case of counter- 
propagating pulses colliding inside the photonic crystal. 
Obviously, it is enough to consider the pulses with iden- 
tical sub-critical intensities and central wavelengths. We 
take the latter to be Ac = 1.064 /im, while the former 
is varied. The results of calculation of the energy excess 
for this case are presented in Fig. |9l It is seen that, at 
low intensities {Am < 2Ao), the excess is negligible which 
means that the pulses propagate independently from one 
another. At higher intensities, a small negative excess 
appears; this means that only a few percent of the en- 
ergy of the pulses is trapped due to their interaction. 
If we increase the pulse amplitude further {Am > 4Ao), 
AW becomes positive, i.e. some additional part of energy 
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second pulse Am2- The amplitude of the first pulse is Ami = 2y4o; the wavelengths Aci = 1.075 /im, Ac2 = 1.064 ^m; the 
time-interval between the pulses At = 5tp; the energy is integrated over the time 300fp. (b) The excess in output energy in 
comparison with the case of single Ara\ and Ami pulses, (c) Profiles of transmitted light for different Aral- (d) Corresponding 
distributions of light intensity inside the photonic crystal at t = SOOtp. 



(less than 20%) is released due to simultaneous presence 
of both pulses inside the structure. 



Thus, there is no evidence of induced trapping in the 
situation of counter-propagating pulses. One can sup- 
pose that the reason is the short interaction time between 
the pulses. Indeed, in the case of co-propagating pulses, 
they travel one after another exchanging energy for a long 
time. On the contrary, the intersection of two subcriti- 
cal counter-propagating pulses is too short-lived to lead 
to any substantial result. The calculations at different 
wavelength show that slowing down of pulses does not 
help to improve this situation. Only at high intensities, 
when the pulses are trapped very soon after launching, 
the traps (which are unstable in this case) seem to be 
sensitive to presence of the second pulse. 



IV. CONCLUSION 

To sum it up, in this paper we have considered the 
interaction of many ultrashort pulses with the photonic 
crystal possessing relaxing cubic nonlinearity. As it is 
known from our previous investigations, there is possi- 
bility of pulse self-trapping in such a structure. First, 
we have analyzed the influence of trapped light on be- 
havior of additional (probe) pulses. It is shown that the 
trap formed by the high-intensity (supercritical) pulse 
can absorb probe pulses with low-intensity, i.e. subcrit- 
ical pulses which do not demonstrate self-trapping on 
their own account. This effect can serve as a peculiar 
absorber for low-intensity pulses. We have also studied 
the changes produced in the trap by incidence of many 
probe pulses. 

The second topic of this paper concerns the induced 
trapping, i.e. the effect of light trapping as a result of 
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collision of two subcritical pulses. We have shown that 
this phenomenon occurs in the case of two co-propagating 
pulses with different velocities: when one of them over- 
takes the other, their interaction leads to effective trap- 
ping of their energy. On the contrary, the coUision of two 
counter-propagating pulses does not allow to observe in- 
duced trapping. 

The results presented in this paper are of general char- 
acter and therefore leave the question of their particular 
realization open. Here we briefly discuss some important 
issues. Wc; have used the simple and well-known Debye 
model of relaxation, though the particular nonlinear me- 
dia can relax according to different laws. This raises an 



interesting question of dependence of the results on the 
relaxation model, thoiigh it is likely that the particular 
model is not as important for the effects calculated as 
the presence of relaxation per se. The choice of model is 
perhaps closely associated with the choice of relaxation 
time. The relaxation times used in our research (up to 10 
fs) make another demand of appropriate materials. Such 
short tni are characteristic, for example, to the electronic 
mechanism and its combinations with other contributions 
to the Kerr nonlinearity. However, our analysis is purely 
phenomenological and does not include these important 
details. 

Another important issue is connected with the value of 
light intensities necessary to obtain the effects reported, 
in particular self-trapping of the pulse. Obviously, inten- 
sities should be very large, putting a question of opti- 
mization of the system in order to reduce the field. Ac- 
cording to this paper, one way to reach this is to use in- 
duced trapping of relatively low-intensity pulses instead 
of one high-intensity pulse. Another possible approach is 
the adjustment of the photonic structure aimed at lower- 
ing the requirements on the materials and pulse param- 
eters. This work is still to be done as well as analysis of 
more general situations. Such situations include taking 
into account of absorption, spatial finiteness of the laser 
beam (we have considered only plane wave approxima- 
tion yet), and other nonlinear contributions (for example, 
light frequency conversion). 



ACKNOWLEDGMENTS 

The work was supported by the Belarusian Foundation 
for Fundamental Research (Grant No. FllM-008). The 
author is grateful to the referees of J. Opt. for interesting 
and inspiring comments. 



[1] S. A. Akhmanov, V. A. Vysloukh, and A. S. Chirkin, Op- 
tics of Femtosecond Laser Pulses (AIP Press, New York, 
1992). 

[2] J. A. Fleck and P. L. Kelley, Appl. Phys. Lett. 15, 313 

(1969). 

[3] V. A. Aleshkevich, S. A. Akhmanov, A. P. Sukhorukov, 
and A. M. Khachatryan, JETP Lett. 13 36 (1971). 

[4] E. G. Hanson, Y. R. Shen, and G. K. L. Wong, Appl. 
Phys.l4, 65 (1977). 

[5] J. A. Fleck and R. L. Carman, Appl. Phys. Lett. 20, 290 
(1972). 

[6] S. Trillo, S. Wabnitz, G. L Stegeman, and E. M. Wright, 

J. Opt. Soc. Am. B6, 889 (1989). 
[7] P. A. Apanasevich, A. A. Afanas'ev, N. B. Gubar, and 

A. L Urbanovich, Laser Phys. 6, 1050 (1996). 
[8] M.-F. Shih, C.-C. Jeng, F.-W. Sheu, and C.-Y. Lin, Phys. 

Rev. Lett. 88, 133902 (2002). 
[9] L Velchev, R. Pattnaik, and J. Toulouse, Phys. Rev. 

Lett.91, 093905 (2003). 
[10] H. Leblond and C. Cambournac, J. Opt. A: Pure Appl. 



Opt. 6, 461 (2004). 
[11] L. Zhang, S. Wen, X. Fu, ,J. Deng, J. Zhang, and D. Fan, 

Opt. Commun.283, 2251 (2010). 
[12] J. M. Soto-Crcspo and E. W. Wright, Appl. Phys. 

Lett.59, 2489 (1991). 
[13] C. Cambournac, H. Maillotte, E. Lantz, J. M. Dudley, 

and M. Chauvet, J. Opt. Soc. Am. B19, 574 (2002). 
[14] S. E. Skipetrov, Opt. Lett.28, 646 (2003). 
[15] X. Liu, J. W. Haus, and S. M. Shahriar, Opt. Com- 

mun.281, 2907 (2008). 
[16] C. Michel, B. Kibler, and A. Picozzi, Phys. Rev. A83, 

023806 (2011). 

[17] A. Armaroli, S. Malaguti, G. Bellanca, S. Trillo, A. de 
Rossi, and S. Combrie, Phys. Rev. A84, 053816 (2011). 

[18] A. C. Ferreira, M. B. C. Costa, A. G. Coelho Jr., C. S. 
Sobrinho, J. L. S. Lima, J. W. M. Mcnezes, M. L. Lyra, 
and A. S. B. Sombra, Opt. Commun.285, 1408 (2012). 

[19] R. A. Vlasov and A. G. Smirnov, Phys. Rev. E61, 5808 
(2000). 

[20] D. V. Novitsky, Phys. Rev. A81, 053814 (2010). 



10 



[21] D. V. Novitsky, Phys. Rev. A84, 053857 (2011). 
[22] D. V. Novitsky, J. Nonlin. Opt. Phys. Mater. 21, 1250010 
(2012). 

[23] D. V. Novitsky, EPL 99, 44001 (2012). 
[24] Z.-M. Meng, F. Qin, and Z.-Y. Li, J. Opt. 14, 065003 
(2012). 



[25] E. D. Palik (cd.). Handbook of Optical Constants of Solids 

(Academic Press, 1998). 
[26] M. J. Weber, Handbook of Optical Materials (ORG Press, 

2003). 



